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ABSTRACT 

We performed a uniform and detailed abundance analysis of 12 refractory elements (Na, Mg, 
Al, Si, Ca, Ti, Cr, Ni, Co, Sc, Mn, and V) for a sample of 257 G- and K-type evolved stars 
from the CORALIE planet search program. To date, only one of these stars is known to har¬ 
bor a planetary companion. We aimed to characterize this large sample of evolved stars in 
terms of chemical abundances and kinematics, thus setting a solid base for further analysis 
of planetary properties around giant stars. This sample, being homogeneously analyzed, can 
be used as a comparison sample for other planet-related studies, as well as for different type 
of studies related to stellar and Galaxy astrophysics. The abundances of the chemical ele¬ 
ments were determined using an LTE abundance analysis relative to the Sun, with the spectral 
synthesis code MOOG and a grid of Kurucz ATLAS9 atmospheres. To separate the Galactic 
stellar populations both a purely kinematical approach and a chemical method were applied. 
We confirm the overabundance of Na in giant stars compared to the field EGK dwarfs. This 
enhancement might have a stellar evolutionary character, but departures from LTE may also 
produce a similar enhancement. Our chemical separation of stellar populations also suggests 
a “gap” in metallicity between the thick-disk and high-a metal-rich stars, as previously ob¬ 
served in dwarfs sample from HARPS. The present sample, as most of the giant star samples, 
also suffers from the B - V colour cut-off, which excludes low-log g stars with high metallic- 
ities, and high-log g star with low-[Ee/H]. Eor future studies of planet occurrence dependence 
on stellar metallicity around these evolved stars we suggest to use a sub-sample of stars in a 
“cut-rectangle” in the logp - [Ee/H] diagram to overcome the aforementioned issue. 

Key words: stars; evolved stars - stars: abundances - stars; planetary systems - techniques; 
spectroscopic - methods: observational 


* Based on observations collected at the Paranal Observatory, ESO (Chile) 
with the Ultra-violet and Visible Echelle Spectrograph (UVES) of the VET, 
under programmes 085.C-0062 and 086.C-0098. 
f E-mail: Vardan.Adibekyan@astro.up.pt 


1 INTRODUCTION 

The precise chemical and kinematic characterization of 
intermediate-mass, evolved stars is very important for differ¬ 
ent fields of both Galactic and stellar astronomy, and the emerging 
field of planetary sciences. 

Many studies observed significant differences in chemical 
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abundances between main sequence dwarf and evolved stars (e.g., 
Friel et aT]|2003| [Jacobson et al.|[2007| jVillanova et al.||2009[ 
Santrich et al.||20f3| l. While these differences for some elements 
might by astrophysical, having a stellar evolutionary character 
(e.g., [Taut vaisiene et al.|[2005| for sodium), several authors how¬ 
ever suggested that the differences may arise also in the analysis, 
being dependent on the particular method and line-list used (e.g., 
[Santos et al.|2009) . Along the same line, one should consider also 
non Local Thermodynamic Equilibrium (non-LTE) effects which 
are stronger for giants than for dwarfs and may have a strong influ¬ 
ence on the analysis (e.g.,[Bergemann et al.|2013[[2014[|Alexeeva[ 
[et al.|2014[ l. 


Understanding the mentioned issues, will not only allow to 
improve of stellar atmosphere models, but also will have very im¬ 
portant implications in several fields of astrophysics. Eor instance, 
it would help us shed light on the statistical and evolutionary prop¬ 
erties of planetary systems around giant stars, e.g., on the possible 
absence of the correlation between stellar metallicity and forma¬ 
tion efficiency of giant planets (e.i 


et al.|[2008 

[Ghezzi et al.[[2010[ [Mortier et al.[[2013| [Maldonado 

etal.|2013 

Jofre et al.l2015r which was found for main sequence 

dwarf stars (e.g.. 

Gonzalez 1997| Santos et al.[2001[ 2004[[Fischer 

& Valentil 

2005 

Johnson et al.[[2010l [Sousa et al.i 201 U [Mortier 


etal.|2013 l. 


Several explanations have been suggested for the aforemen¬ 
tioned lack of metallicity enhancement for giant stars hosting a gi¬ 
ant planet. Higher stellar mass of giants may compensate the lack 
of metals (e.g.,[Ghezzi et al.|2010|; possible spectroscopic analysis 


issues in giant stars (e.g., Hekker & Melendez|2007[ Santos et al.[ 
[2009[l; selection biases in giant star samples iMortier et al.|2013|>. 


However, one should note that some studies reported an enhanced 
metallicity of giant stars with planets, but with small samples of 


planet hosts JHekker & Melendez 

2007 

IQuirrenbach et al.|2011|l. 

We refer the reader tolAlws et al.j 

2015 

and reference therein) for 


more detailed review on the topic. 

In this paper we focus on the chemical and kinematic proper¬ 
ties of a sample of 257 field giant stars which are observed within 
the context of the CORALIE extrasolar planet search program. The 
main characteristics of the sample along with the homogeneously 
derived stellar atmospheric parameters are presented in a parallel 
paper ( [Alves et al.|2015^ . The uniform chemical analysis of these 
giant stars is very important to explore the specific chemical re¬ 
quirements for the formation and evolution of planetary systems 
around them. The paper is organized as follows: in Sect. 2, we 
briefly introduce the sample used in this work. The method of the 
chemical abundance determination and analysis will be explained 
in Sect. 3. The distinction of different Galactic stellar populations 
and kinematic properties of the stars are presented in Sect. 4. Then, 
after discussing the metallicity distribution of the stars in Sect. 5, 
we summarize our main results in Sect. 6. 


2 SAMPLE DESCRIPTION AND STELLAR 
PARAMETERS 

Our sample comprises of 257 G- and K-type evolved stars that are 
being surveyed for planets in the context of the CORALIE (|Udry| 
[et al.|2000l( extrasolar planet search program. High-resolution and 


Table 1. The slope, correlation coefficient, and the significance of the 
[X/Fe] linear trends with the Teff. 


Elem 

slope 


R2 

N 

z-score 

Nal 

1.62±0.36 X 

10-4 

0.72 X 10-4 

256 

4.3 

Mgl 

-2.51±0.38 X 

10-4 

0.14 X 10“ 

257 

6.1 

All 

-2.69±0.29 X 

10-4 

0.24 X 10“ 

257 

7.9 

Sil 

-2.39±0.16 X 

10-4 

0.45 X 10“ 

257 

11.0 

Cal 

0.61±1.80 X 

10-® 

0.44 X 10-5 

257 

0.3 

Scl 

-3.00±0.34 X 

10-4 

0.23 X 10“ 

255 

7.6 

Sell 

-0.78±0.23 X 

10-4 

0.43 X 10-4 

257 

3.2 

Til 

0.81it2.94 X 

10-5 

0.29 X 10-5 

257 

0.2 

Till 

-0.44±0.34 X 

10-4 

0.64 X 10-5 

257 

1.2 

VI 

-1.79±0.29 X 

10-4 

0.13 X 10“ 

256 

5.9 

CrI 

0.56±0.10 X 

10-4 

0.98 X 10-4 

257 

4.9 

CrII 

1.33±0.34 X 

10-4 

0.57 X 10-4 

246 

3.6 

Mnl 

1.33±0.34 X 

10-4 

0.57 X 10-4 

247 

2.3 

Col 

-0.32±0.03 X 

10-5 

0.33 X 10“ 

257 

9.2 

Nil 

-0.78±0.11 X 

10-4 

0.16 X 10“ 

257 

6.3 


high signal-to-noise (S/N) spectra were obtained using the UVES 
spectrograph. Precise stellar parameters for the entire sample were 
determined in [Alves et al.[ ( [20T5) ( by using the same spectra as we 
did for this study. The spectroscopic stellar parameters and metal- 
licities were derived by imposing excitation and ionization equi¬ 
librium. The spectroscopic analysis was completed assuming LTE 
with a grid of Kurucz atmosphere models ( [Kurucz|1993t, and th e 
2002 version of the MOOCj^adiative transfer code i Sneden 1973 i. 
We refer the reader to [Alves et al.[ ( [20T5l l and |Sousa ( [2014^ for de¬ 
tails. 


[Alves et aL| (|2015|( derived the atmospheric parameters by us¬ 
ing three different line-lists of Eel and Eell ( [Hekker & Melendez[ 
[2007[[Sousa et al.|2008[[Tsantaki et al.|2013[ l. Whilst showing that 

the use of different line-lists gives compatible results, the parame¬ 
ters derived following [Tsantaki et al.[ ( [2013[ were adopted, so we 
also do for the rest of the present paper. 

The stars in the sample have effective temperatures 4700 < 
Teff < 5600 K, surface gravities 2.2 < \ogg < 3.7 dex, microtur¬ 
bulence 1 3.2 km s“^and they lie in the metallicity range 

of -0.75 < [Fe/IU < 0.3 dex. 


3 CHEMICAL ABUNDANCES 

For the abundance derivation we closely followed the method de¬ 
scribed in [Adibekyan et alT| ( |2012| (. 

3.1 Selection of the lines and abundance derivation 

The initial line-list and the atomic data were taken from Adibekyary 
[et al.[pOT2) and [Neves et al.[ ( [200^ . [Neves et al.[ ( |2009 1 provieded 
the astrophysical (callibrated) oscillator strength and solar equiv¬ 
alent widths of the lines. Since the spectra of cool evolved stars 
are more line crowded (which cause strong blending) compared to 
their unevolved hotter counterparts, we aimed to carefully select a 
subset of unblended lines from [Adibekyan et al.[ ( |20f2l l. For this 
purpose as a reference we used a very high S/N and high resolu¬ 
tion archival spectrum of the K-type giant Arcturus observed with 


^ Indeed, [Reffert et al.[ [2014} claims a strong evidence for a planet- 
metallicity correlation for giant planet host stars. 


^ The source code of MOOG can be downloaded at 
http://www.as.utexas . edu/~chris/moog.html 
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Figure 1. [Crl/Crll], [ScI/ScII], and [Til/Till] as a function of atmospheric parameters for our sample of evolved stars (black points) and for the sample of 
FGK dwarf stars from |Adibekyan et al.H2Ql^ (grey dots). 


the NARVAL spectrograph ( [Mortier et al.|2013) l. We measured the 
equivalent widths (EWs) of the selected lines both manually, using 
a Gaussian fitting procedure within the I RAIp]splot task , and 
automatically, by using the AREI^code i Sousa et al.|2007 1 . We 
calculated the mean relative difference ((EWaabs - EW/_raf)/ 
EW/baf) and standard deviation of the relative difference of the 
EW measurements and applied 2(T-clipping. We repeated this pro¬ 
cedure a second time after the outliers were excluded. Einally, 118 
lines out of 164 were left that show a relative difference in EW of 
less than 15%. These lines were once again checked by eye within 
IRAF to make sure that they are not blended and hence the corre¬ 
spondence between the EW measurements is not by chanc^ 

After selecting the isolated lines, the abundances for 12 ele¬ 
ments (Na, Mg, Al, Si, Ca, Ti, Cr, Ni, Co, Sc, Mn, and V) were 
determined using a LTE analysis relative to the Sun with the 2010 
version of the MOOG ( |Sneden|1973^ and a grid of Kurucz ATLAS9 
plane-parallel model atmospheres with no a-enhancement. The ref¬ 
erence abundances used in the abundance analysis were taken from 
[Anders & Greves^ ( |1989^ . We note that our analysis is differential 
and the source of the reference abundances is not crucial. For the 
automatic EW measurements we used ARES for which the input 
parameters were the same as in jSousa et al.| (|2008|l and the rejl pa¬ 
rameter is calculated following the procedure discussed in [Mortier] 
[et al.[ ( |20T3l l. The EWs of all the lines used in the derivation of the 
abundances for all studied stars is available at CDS. 

The final abundance for each star and element was calculated 
to be the average value of the abundances given by all lines detected 
in a given star and element. Individual lines for a given star and 
element with a line dispersion more than a factor of two higher than 
the rms were excluded. In this way we avoided the errors caused 
by bad pixels, cosmic rays, or other unknown effects. A sample of 


® IRAF is distributed by National Optical Astronomy Observatories, op¬ 
erated by the Association of Universities for Research in Astronomy, Inc., 
under contract with the National Science Foundation, USA. 

^ The ARES code can be downloaded at 
http://www.astro.up.pt/'~sousasag/ares 
The line-list is available at the CDS 


our results for three stars is presented in Table and the complete 
results are available at the CDS. 


3.2 Uncertainties 

Since the abundances were determined via the measurement of 
EWs and using already determined stellar parameters, the errors 
might come from the EW measurements, from the errors in the 
atomic parameters, and from the uncertainties of the atmospheric 
parameters that were used to make an atmosphere model. In ad¬ 
dition to the above-mentioned errors, one should add systematic 
errors that can occur due to NLTE or granulation (3D) effects. To 
minimize the errors, it is very important to use high-quality data 
and as many lines as possible for each element. 

We followed |Adibekyan et al.|(2012[ l for the calculation of the 
errors. In short, we first varied the model parameters by an amount 
of their one-sigma errors available for each star and calculated the 
abundance differences between the values obtained with and with¬ 
out varying the parameter. Then we evaluated the errors in the abun¬ 
dances of all elements [X/H], adding quadratically the line-to-line 
scatter errors and errors induced by uncertainties in the model at¬ 
mosphere parameters. In cases when only one line used to derive 
the abundances a typical 0.1 dex error for line-to-line scatter was 
assumed. For our sample stars the errors induced by uncertainties 
in the parameters of model atmosphere varies from about 0.02 dex 
(for Sil) to « 0.06 dex (for VI) and in general are smaller than the 
line-to-line scatter errors. The final errors for the studied elements 
are smallest for All (« 0.04 dex) and largest for VI (~ 0.14dex). 


3.3 Testing the validity of the stellar parameters 

Although [Alves et al.[ ( [2015) have shown that the stellar parameters 
in general agree very well with the literature data, the consistency 
does not always imply correctness. Moreover the stellar parame¬ 
ters were derived by completing an LTE abundance analysis and 
by using only iron lines. To check the validity limit of the adopted 
methodology in terms of stellar parameter ranges we tested our re¬ 
sults in two ways (see also[Adibekyan et al.|2012^. 
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Figure 2. [X/Fe] vs. Teff plots. The black dots represent the stars of the sample and the gray small dots represent stars from |Adibekya n et al.H201^ . The blue 
and red solid lines depict the linear fits of the current data and the data from |Adibekyan et al.H20l'^ , respectively. The blue dashed line is the fit of our data 
after correcting for the trend with Teff. Each element is identified in the upper right corner of the respective plot. 


Table 2. Sample table of the derived abundances of the elements, rms, total error, and number of measured lines for each star. 


Star 

... [Sil/H] 

rms 

err 

[SiI/H]*o^^ 

n 

[Cal/H] 

rms 

err 

n 

HD47001 

-0.20 

0.09 

0.09 

-0.26 

14 

-0.25 

0.05 

0.06 

11 

HD73898 

-0.30 

0.03 

0.03 

-0.28 

14 

-0.32 

0.04 

0.05 

11 

HD16815 

-0.16 

0.07 

0.07 

-0.20 

15 

-0.25 

0.04 

0.06 

12 ... 


Notes: * The [X/H] abundances after correction for the T^ff trends. 
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Figure 3. [X/Fe] vs. [Fe/H] plots. The black dots represent the stars of the sample and the gray small dots represent stars from |Adibekyan et al.H201^ with 
Teft = Tq ±500 K. The red circle and blue square show the average [X/Fe] value of stars with [Fe/H] = 0.0±0.1 dex. Each element is identified in the upper 
right corner of the respective plot. 


First, in Fig.[T]we plot the [CrI/CrII], [ScI/ScII], and [Til/Till] 
as a function of the stellar parameters to ensure that the ionization 
equilibrium enforced on the Fell lines jAlves et al.|2015^ is accept¬ 
able to other elements. For comparison, the field FGK dwarf stars 
from |Adibekyan et al.| ( |201^ are presented. Most of the trends are 
nearly flat around zero in contradiction to their unevolved coun¬ 
terparts for which a gradual increase with decreasing TeS was ob¬ 
served. For our stars only an increase of [ScI/ScII] ratio can be seen 
with the decrease of Tes. However, the results obtained for Scl and 
CrII should be considered with caution since their abundances have 
been derived by using only one line. From the figure one can notice 
a small offset from zero for [Til/Till] ratio and [ScI/ScII]. These 


positive offsets probably do not have relation to the NLTE effects 
as discussed in |Bergemann| ( [20TT) for [Til/Till] and still need to be 
understood. 

ITsantaki et al.| ( |2013) showed that by correcting stellar pa¬ 
rameters (mainly Teft). using carefully selected line-list especially 
designed for cool stars, the observed trends of [XI/XII] with stellar 
parameters get flatter. For example, an overestimation of Tes for 
cool stars might cause of the trends. The dependence of the abun¬ 
dances of ionized and neutral species on the surface gravity is also 
discussed in |Mortier et al.| ( (2013[ l. 

Another way of testing the stellar parameters is to plot [X/Fe] 
against TeS (Fig. |^. For the comparison the dwarf sample of 









































6 V. Zh. Adibekyan et al. 


0.4 

0.3 

0.2 

0.1 

0 


- 0.1 


• 0.2 

0.4 

0.3 

0.2 


0 



- 0.1 


• 0.2 

0.4 

0.3 

0.2 

0.1 

0 


- 0.1 


- 0.2 


' 1 ' 1 ' 1 ' 1 ' 1 

' 1 ' 1 ' 1 ' 1 ' 1 

' 1 ' 1 ' 1 ' 1 ' 1 

•• .. 

./ ^ V ..•••• • 

• 

• • 

.. . . \ - 

, •••• • 

i - 


_ X=5645.61 A; Xi=4.93 eV_ 

1 1 1 1 1 1 1 1 1 1 

_ X=5701.ll A; Xi=4.93 eV_ 

1 1 1 1 1 1 1 1 1 1 

_ X=5684.49 A; x,-4.95*eV_ 

1 1 1 1 1 1 1 1 1 1 

1 1 1 1 1 1 1 1 1 1 

• 

r. . 

• W 

• *n5^ 

1 1 1 1 1 1 1 1 1 1 

• 

L*t ** * 

.1 1 1 1 1 1 1 1 1 1 

- • • • - 
•*•* • 

• ^ •- 

•• 

_ X = 6145.02 A; Xi=5 62 eV_ 

1 1 1 1 1 1 1 1 1 1 

_ X=6142.49 A; Xi = 5-62 eV_ 

il 1 1 1 1 1 1 1 1 1 

_ X = 5753.64 A; ;^,=5.62 eV_ 

1 1 1 1 1 1 1 1 1 1 

1 ' 1 .L 1 ' 1J 1 

^ - • 

• % • 

• * •fflPsr 

** • K. «. 

•' 1 ' J 1 ' 1.' *1 

• •• • •- 

V - 

1 1 1 1 1 1 1 1 1 1 

% 

• V • . - 

J_ ^ 



_ \ = 6195.46 A; = eV_ 

. 1 . 1 . 1 . 1 . 1 

_ X=6527.21 A; Xi = 5.B7 eV_ 
. 1 . 1 . 1 . 1 . 1 

_ X = 6741.63 A; y, = 5.98 eV_ 

.. 


-0.6 -0.4 -0.2 0 0.2 -0.6 -0.4 -0.2 0 0.2 -0.6 -0.4 -0.2 0 0.2 

[Fe/H] 


Figure 4. [Si/Fe] vs. [Fe/H] plots for Si lines of different excitation potentials. The black dots represent the stars of the sample and the gray small dots represent 
stars from [Adibekyan et'n!]j2012) with Tgff = Tq ±500 K. The wavelength of each line and excitation energy of the lower energy level (xi) is identified in 
the lower left comer of the respective plot. 


|Adibekyan et al.|(2012| l is also presented. Stellar evolutionary mod¬ 
els do not suggest significant trends of these ratios with Teff. How¬ 
ever, for several elements we detected significant trends. To eval¬ 
uate the significance of the trends we performed a linear fit and 
followed the procedure described in |Figueira et al.| ( [2013^ . In short, 
first, we obtained the zero-centered distribution of the correlation 
coefficient by randomly bootstrapping (building random samples 
by shuffling the parameters among the observed set of parameters) 
the observed data pairs 10"^ times. Then we calculated the corre¬ 
lation coefficient for each of these uncorrelated data sets and then 
the average and standard deviation of these values. By assuming 
a Gaussian distribution for R (correlation coefficient), we calcu¬ 
lated the probability that the R of the original dataset was obtained 
by pure chance. The significance of the trends and the slopes are 
presented in Table [T] From the figure, one can see that for most 
cases the trends are less steep compared to those observed for the 
dwarf]^ which in turn speaks about the correctness of the stellar 
parameters used to derive the abundances. 


° Note that the stellai' parameters for the dwarfs were not derived by using 
the |Tsantaki et al.| j2013) line-list which is especially designed for cool 
stars. 


[Adibekyan et al.|f2012^ already discussed several possible rea¬ 
sons for the observed trends of [XFXII] with stellar parameters, 
and [X/Fe] with TcaQ and concluded that the observed trends are 
probably not an effect of stellar evolution, and uncertainties in at¬ 
mospheric models are the dominant effect in the measurements. 
The authors afterwards removed the Ten trend as it was done also 
in other works (e.g., [Valenti & Fischer||2005| [Petigura & Marcy[ 

|20TT] |. 

Since by fitting the data and simply subtracting the fit would 
force the mean [X/Fe] to zero (which is an nonphysical situation), 
[Adibekyan et al.[ \1Q\2) added a constant term chosen so that the 
correction is zero at solar temperature. In our case the stars are 
cooler and their temperatures do not reach the solar temperature so 
we decided to apply another approach. In this case, the constant 
term was chosen so that the correction is zero at = 4960 K, 
which is the mean temperature of our sample stars. However, we 
appreciate the fact that this approach and the choice of the con¬ 
stant term is arbitrary. For this reason we decide to use the original 
(before detrending) chemical abundances for the rest of our study. 

^ Note that the [Adibekyan et al.j ^2012^ sample essentially consists of 
dwarfs. 
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Figure 5. Toomre diagram for the entire sample. The left and right panels show the separation of the stellar groups according to the|Bensby et al.|^2003| B03) 
and|Robin et al.|j2003| R03) prescription, respectively. The symbols are explained in the figure. 


The dependence of [X/Fe] on stellar gravity and microturbu¬ 
lence and metallicity is shown in Fig. |Bl[ Fig. |B2[ and Fig.[^ re¬ 
spectively. For most of the species we did not observe a trend with 
and log g, and some of the observed trends probably reflect the 
correlation of TeS with other stellar parameters (see Fig. 1^. 

As a final check, we compare our derived abundances with 
those obtained by |Liu et al.| ( |2007) . We note that this is the only lit¬ 
erature source where we find enough stars (14 stars) in common to 
compare. We found very good agreement for all the species except 
vanadium: /S.[Na/H] = 0.02±0.12, A[Mp/iT] = -0.06±0.12, 
A\Al/H] = 0.05 ± 0.04, A[Ca/H] = 0.03 ± 0.06, /S.[Si/H] = 
0.03±0.03, A[rf/T/] = 0.04±0.08, A[Af/H] = -0.01±0.05, 
and A[l//F/] = 0.19 ±0.050 

The [X/H] abundances of all the stars before and after correc¬ 
tion (if the significance of the correlation is above 3 (t) for the Teff 
trends are available at the CDS (see also Table]^. 


3.4 [X/Fe] vs. [Fe/H] relation. Galactic chemical evolution 


The [X/Fe] vs. [Fe/H] relation plot is traditionally used to study the 
Galactic chemical evolution because iron is a good chronological 
indicator of nucleosynthesis. In Fig. 0 we present the dependence 
of [X/Fe] on met allicity for our sample of g iant stars and for FGK 

In the figure we also 


Adibekyan et al. 


m ple ot 

1 (2012!- 


dwarf stars from 
showed the average value of [X/Fe] for stars in the metallicity range 
of 0.0 ±0.1 dex, where the Galactic chemical evolution effects are 
small. As one can see, for all the elements the general behavior of 
[X/Fe] with the metallicity is similar for giant and dwarf stars, and 
reflects the Galactic chemical evolution in the solar neighborhood. 
However, one can also clearly notice that, for some elements (Co, 
Na, V, Mn, Al, and Si) while the Galactic chemical evolution trends 
are similar, they are shifted: for giant stars having higher [X/Fe] 
values at a fixed metallicity. The largest offset is seen for Na and 


» A[X///] = [X/H]our - [X/H]theirs- 

^ Only stars with Teff = Tq ±500 K are presented, because of the highest 
accuracy in the parameters and chemical abundances in these stars. 


Mn, and a bit less in Si and Al. In general, Na and Al are not good 
tracers of chemical evolution and affected by internal mixing pro¬ 
cesses in the stars. The Mn abundance was obtained by using only 
one line and it should be considered with caution. Moreover, the 
scatter in [Mn/Fe] is very high, indicating unrealistic abundances 
for some fraction of the stars. 

Overabundances of sodium and aluminum (also silicon in 
some cases) in open cluster giants (compared to the abundances 
of dwarfs) were already observed by several authors (e.g., Friel| 
let al.|200^|2005]|Tautvaisiene et aLl|2005| [Jacobson et al.|2007| 

IVillanova et al.|2009[|Santrich et al.|2013| >. In most of these stud¬ 
ies, the trends were explained as a stellar evolutionary effect, due 
to the deep mixing produced by the hydrogen burning cycle, af¬ 
ter stars have left the main sequence. For a complete picture, one 
should perform thorough analysis taking into account the non-LTE 
effects which are stronger for giants stars and also the systematic 
errors which might arise due to particular spectroscopic analysis 
method used. For example, it is well known that sodium lines suf¬ 
fer from non-LTE effects which lead to an overestimation of the 
Na abundances (e.g., |Alexeeva et al.|2014] . In our analysis we used 
two sodium lines (at 6154.23 Aand 6160.75A) which were studied 
for non-LTE effects in | Alexeeva et al.H2014| l. The average EWs of 
these lines were ~ 70 mA for 6154.23 A, and ~ 80 mA for the 
6 I 6 O. 75 A line. According to Alexeeva et al., the non-LTE correc¬ 
tion for our stars should be from -0.1 to -0.15 dex, which is close 
to the difference in [Na/Fe] between giants and dwarfs observed in 
this study. 

The difference in Al abundances ([Al/Fe]) between giants and 
dwarfs obtained for solar-metallicity stars is not large (0.07 dex 
- about la scatter), but seems to increase at lower metallicities. 
However it should be noted that direct comparison of the abundance 
ratios at lower metallicities is not straightforward, since the Galacic 
chemical evolution effects and the relative fraction of thin and thick 
disk stars can be dominant. Several authors studied the non-LTE 
effects on the formation of Al lines (e.g. |Baumueller & Gehrenj 
|1996|p'997]|Menzhevitski et al.|2012] . They showed that the non- 
LTE correction of the Al abundances, derived from the subordinate 
doublet AA 6696.03, 6698.68 A, is very small at solar metallicities. 
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does not depend strongly on the surface gravity and only significant 
at temperatures above 6500 K ( |Menzhevitski et al.|2012| ). 

The next element for which we obtained small, but system¬ 
atic difference between giants and dwarfs is Si. The abundance 
of Si is not expected to be affected be extra mixing processes in 
the stars. The few studies of the Si abundances taking into account 
the non-LTE deviations showed that the effect is significant only 
at very low metallicities (e.g. |Bergemann et al.|2013) and the non- 
LTE correction for Si of the Sun is about -0.05 dex |Sukhorukov ^ 
|Shchukina|2012^ . Since in this study for the Al abundance deriva¬ 
tion we used several lines with different excitation potentials (EP), 
which means different atmospheric layers of formation and hence 
different sensitivities to non-LTE deviations, we decided to ana¬ 
lyze [Si/Ee] against [Fe/H] for each individual line. In Fig. we 
plot [Si/Fe] against [Fe/H] for 9 Si lines with the lowest, interme¬ 
diate and highest excitation energy of the lower energy level (xi)- 
At the first glance it looks like the lines with the highest xi show 
the highest deviations from the abundances derived for the FGK 
dwarfs. However, as can be seen in the middle panel of the plot the 
three lines with exactly the same xi show different behavior, for A 
5753.64 A showing the largest difference. 

Fig. El shows that the picture is complex and probably sev¬ 
eral process (NLTE, unresolved blends) are acting and affecting 
the abundances at the same time. To select the “besf’ lines i.e. 
lines which give similar average abundances to that obtained for 
the dwarfs, for all the lines we calculated the average [Si/Fe] abun¬ 
dance ratio obtained for all the giants with solar-metallicity ±0.1 
dex and compared that with the average [Si/Fe] obtained for dwarfs 
( [Adibekyan et al.|2012| l with meallicities in the range of 0.0 ±0.1 
dex. Then we used the rms (scatter) of the [Si/Fe] calculated for the 
dwarf j^to quantify the observed differences (in nxo-). We found 
five (out of 15) Si lines which give [Si/Fe] abundance similar to that 
obtained for the dwarfs (less than 2a difference). The [Si/Fe] abun¬ 
dance obtained by averaging the abundances of the mentioned five 
lines against [Fe/H] is presented in the Fig. |B3| We note that we do 
not claim that these selected lines are not affected by non-LTE ef¬ 
fects of unresolved blends, but they provide abundances similar to 
that obtained for dwarfs, which probably means that they are more 
realistic. 

We repeated the aforementioned procedure for all the lines for 
each element and calculated the difference in [Xjine/Fe] between 
giants and dwarfs for each line. This differences,in nxcr, is pre¬ 
sented in the last column of the line-list table (available at CDS). 
The [X/Fe] vs. [Fe/H] relation obtained by using only the “best” 
lines (less than 2cr difference) is shown in Fig. |B3| The correspond¬ 
ing [X/H] abundances are available at CDS. 


4 KINEMATICS AND STELLAR POPULATIONS 


It is becoming increasingly clear that a separation of the Galac¬ 
tic stellar components based only on stellar abundances is supe- 


rior to kinematic separation (e.g.„ Navarro et al.|2011| Lee et al. 

2011 

Adibekyan et al.[2011|[Liu & van de Ven[2012 Recio-Blanco 

et al. 

20 14^, because chemistry is a relatively more stable prop- 


erty of a star than its spatial positions and kinematics. However, as 
mentioned above, some changes in abundances of some elements 


are expected when the stars are evolving and leaving the main se¬ 
quence. In this analysis, to separate the thin- and thick disk stellar 
components, we used the position of the stars in the [a/Fe] - [Fe/H] 
plane (here a refers to the average abundance of Mg, Si, and Ti), 
but separately also a kinematics approach is applied. 

The space velocity components for 183 stars out of 257 were 
derived with respect to the local standard of rest (LSR), adopting 
the standard solar motio n (Ur , Vq, W©) = (11.1, 12.24, 7.25) km 


of 


Schdnrich et al. 


1 2010 1 . For the remaining 73 stars we did 


not calculate the velocities because of the deficit of astrometric lit¬ 
erature data. The radial velocities, parallaxes and proper motions 
were taken from the SIMBAD Astronomical Database I Combin¬ 
ing the measurement errors in the parallaxes, proper motions, and 
radial velocities, the resulting average errors in the U, V, and W 
velocities are of about 2-3 km s“^. 

To assess the likelihood of the stars being a member of dif¬ 
ferent Galactic populations, we followed [Reddy et al.| ( |2006| l. The 
probabilities that the stars belong to different stellar populations 
were calculated, having adopted both the [Bensby et al.[ ( [2003 |l and 
Robin et al. (2003^ population fractions. We refer the reader to 


Adibekyan et al. 1 2012[( for the details of the computation. The 


Galactic space velocity components and the probabilities to assign 
the stellar population to which the stars belong are available at the 
CDS. 

According to the [Bensby et al.[j2003^ criteria, among the 183 
stars, we have 176 (96%) stars from the thin disk, 5 from the thick 
disk, and 2 stars are considered to be transition stars that do not 
belong to any group. Adopting the criteria from [Robin et al.[j2003[ ( 
gives 177 (97%) thin disk stars, 5 star with kinematics suggest¬ 
ing a thick/thin disk transition, and one star with a classification of 
thick-disk/halo transition object. The distribution of the stars in the 
Tommre diagram is shown in Fig.[^ 

As mentioned above, in addition to the difference in their kine¬ 
matics, the thin- and thick disk stars are also different in their a 
content at a given metallicity (e.g., [Fuhrmann[1998[[2008] (. This di¬ 
chotomy in the chemical evolution allows one to separate different 
stellar populations. 

The [a/Fe] versus [Fe/H] plot for the sample stars along with 
the dwarf stars from jAdibekyan et al.[j201^ with Teff = r0±5OO 


K is depicted in Fig. f ^ As one can see from the figure the two 


samples show similar trends, with giant stars having on average 
higher [a/Fe] values at a fixed (low) [Fe/H]. This observed differ¬ 
ence might arise from our assumption of LTE line formation. The 
non-LTE effects are stronger for metal-poor stars, but these effects 
depends also on other stellar parameters (e.g. gravity and temper¬ 
ature) and also they different for different elements and they differ 
from line to line. Thus, for fully understanding of the main reason 
of the observed abundance difference between giants and dwarfs a 
complete non-LTE analysis is needed. 

Our chemical separation of the Galactic disks suggest that 
23 stars (9%) in the sample show enhanced a abundances. In 
[Adibekyan et aL[ j2011^ and [Adibekyan et al.[ ( [2013[ l the high-a 
stars were separated into two groups with a gap in both [a/Fe] and 
metallicity. It is interesting to see that the gap in [Fe/H] for high-a 
stars can be also seen in our sample at the same metallicity (« -0.2 / 
-0.3 dex). Following the same logic and definitions as injAdibekyanj 
leTaLH^OT^ , the 10 stars with enhanced [a/Fe] and [Fe/H] below 
-0.3 dex can be classified as thick-disk stars, and the remaining 13 


The scatter obtained for the dwaifs by averaging the abundance of many 
lines is more realistic than the scatter obtained from one spectral line for 
giant stars. 


http://simbad.u-strasbg.fr/simbad/ 

The chemical dissection of the disks is presented in the appendix. 
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Figure 6. [o/Fe] vs. [Fe/H] for the current sample (black dots) and for the 
stars from [Adibekyan et al.|j20l^ with = FQitSOO K (gray small 
dots). The separation between the thick- and thin-disk stars for the two sam¬ 
ples are presented in black and gray dashed lines. 


stars as high-a metal-rich stars (hamr). With this definition we see 
that 4% of the stars belong to the Galactic thick disk, as the kine¬ 
matic separation was suggesting. 

We note that the current sample is small and we will avoid of 
a definite conclusion about the existence of the mentioned “gap" 
at [Fe/H] ~ -0.3 dex and the distinction of the two a-enhanced 
metal-poor and metal-rich populations. However, the fact that the 
two different homogeneously analyzed samples (the current one 
and the one from [Adibekyan et al.| ( |20lT| l) show quite similar fea¬ 
tures probably is more than just a hint about the existence of the 
hamr stars as a distinct stellar family. Moreover, recent study of 
an inner disk metal-rich open cluster, Berkeley 81, shows that the 
stars are enhanced in a element Magrini et al. (2014, in prep.), 
thus confirming that the hamr stars have inner disk origin as sug¬ 
gested by [Adibekyan et aT] ( |2013|). However, we want to note that 
no similar ”gap“ was found in Bensby et al.[ ( [20l'4l > where the au¬ 
thors suggested that the hamr stars present the metal-rich tail of 
the thick disk. As mentioned in [Bensby et al.[p014^ , a large sam¬ 
ple with well-controlled selection function (e.g., Gaia-ESO survey 
- [Gilmore et al.[ ( [?012[ l) would help us to understand the real nature 
of the hamr stars. 


5 METALLICITY DISTRIBUTION 

As mentioned above, several authors tried to understand the reason 
why the apparent giant-planet - metallicity correlation does not ex¬ 
ist for evolved stars. As recently suggested by [Mortier et al.[ ( [2013[ l, 
a possible reason might be a selection bias due to S - V colour 
cut-off. 

In Fig. 1^ we plotted the relation between stellar metallicity 
and surface gravity. For the comparison, the dwarf stars sample 
from [Adibekyan et al.[j2012) is also presented. From the figure one 
can easily see that the giant stars sample lacks high-metallicity and 
low-gravity stars, and also low-metallicity and high-gravity stars. 
This is again probably because of the selection criteria used to de¬ 
fine the sample. 

To avoid the issues related to the selection effects, an unbiased 
giant sample with no colour cut-off and homogeneously derived 
parameters is needed that is systematically searched for planetary 


companions. However, it is still possible to overcome the effect of 
the B - V colour cut-off if one considers, for example, only stars 
in the "cut rectangle" shown in Fig. (red rectangle), where the 
stars are equally distributed. However, these "cut rectangles" will 
consist of stars with narrower ranges of metallicities (from -0.25 to 
0.15 dex in the example of Fig. [^, which is also an issue since the 
giant-planet - metallicity correlation is more pronounced at high 
metallicities (at least for dwarf stars). 

In the right panel of Fig. |7] we show the metallicity distribu¬ 
tion of giant and dwarf stars where narrower [Fe/H] distribution of 
giants is apparent. The figure also shows that the two distributions 
are peaked at similar metallicities, close to the solar value. The me¬ 
dian (and its standard deviation) of the metallicity distributions of 
giant and dwarf stars are -0.05 (0.18) dex and -0.10 (0.33) dex, re- 
spectivel}^^ Several studies have already observed this tendency 
of evolved stars lacking the metal-rich and very metal-poor tails 
(e.g., [Taylor & Croxall[[2005| [Takeda et al.[[2008| [Luck & Heiter[ 

|2007l[Ghezzi et al.|2010t. 


The stars in this sample have stellar masses between 1.5 and 


4.0 Mq <Alves et al.[[2015 

1 , and hence should be on average 

younger than the dwarfs from 

Adibekyan et al.[([2012J. The younger 

age together with the age - metallicity dispersion relation (e.g., [da 

Silva et al.|2006[[Haywood|2008|[Casagrande et al.|2011|[Maldon- 

ado et al.|2013| might explai 

n the narrower [Fe/H] distribution of 

)stly local since they do not have time 
^Wang & Zhao|2013[[Minchev et al. 
e disk is expected to make the metal- 
loes not change the mean abundance 
see in Fig.|^ This is because mostly 

the giants. Young stars are me 
to migrate within the Galaxy 
2013[l. Radial migration in th 
licity distribution wider, but 
([Wang & Zhao|2013[l, as we 


massive stars contribute to the chemical enrichment of the inter¬ 
stellar medium and they contribute mainly around their birth places 
because of their very short lifetime. The lack of very metal-rich gi¬ 
ants can be understood along the same migration process, most of 
the old stars which migrate would come from the inner, metal-rich 
disk ( [Wang & Zhao[20T^[Minchev et al.[2013[ >. 

In addition to the aforementioned astrophysical explanation, 
we would like to note again the selection effects which may arise 
in evolved star samples due to B -V colour cut-off. This selection 
bias may also make the metallicity distribution narrower. 


6 SUMMARY OF THE RESULTS 

We have carried out a uniform abundance analysis for 12 refractory 
elements for a sample of 257 field G-, K-type evolved stars that 
are being surveyed for planets using precise radial-velocity mea¬ 
surements with the CORALIE spectrograph. The abundances were 
derived using a carefully selected line-list and are based on the pre¬ 
cise spectroscopic parameters derived by [Alves et al.[ ( (^15[ l using 
the same spectra as were used in the present study. 

We found that for all the elements Galactic chemical evolu¬ 
tion trends are similar for giant and dwarf stars, while for some 
species [X/Fe] values are shifted towards higher values at a fixed 
metallicity. Our LTE analyis confirms the overabundance of Na in 
giant stars compared to the field FGK dwarf stars from [Adibekyan[ 
[et al.[ ( [^12[ (. This overabundance may have a stellar evolutionary 
character, even though the possible departures from non-LTE may 
produce an enhancement of a similar degree ([Alexeeva et al.[20T4|(. 


We note that the standard deviation of the median is calculated as 
1.25 *(t, where a is the standard deviation of the distribution. 
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Figure 7. Left panel: [Fe/H] vs. log g for the current sample (black dots) and 
for the stars from | Adibekyan et al.lpOlT) (gray dots). The two black dashed 
lines were drawn by eye and show the bia.ses in the samples due to the B - V 
cut-off. Right panel: The metallicity distribution of the two aforementioned 
samples. The distribution of the giants stars (gray line) was multiplied by 2 
for the better visual comparison. The median and its standard deviation is 
also presented for metallicity distributions of both giants and dwarfs. 


We showed that an observed small overabundance of Si compared 
to the field FGK dwarf vanishes when a shorter, carefully select 
line-list is used. 

To separate Galactic stellar populations, we applied both a 
purely kinematical approach and chemical method. Our chemical 
separation suggests that 91% of the stars, being a-poor, belong to 
the thin disk and the remaining 9% of the stars show enhanced a- 
element abundances at a fixed [Fe/H]. This sample (while being not 
very large) also suggests a ”gap“ in [Fe/H] for high-a stars as ob¬ 
served in [Adibekyan et ahj l lTOll^ . Following the definition of the 
last authors 4% of the stars were classified as thick-disk members 
(being metal-poor) and 5% as hamr stars. 

The metallicity distribution of the giant stars is shown to be 
narrower than that of their non-evolved dwarf counterparts (see also 
[Taylor & Croxall|2005[ [Takeda et al.|2008) , but peaked at almost 
solar metallicity as in case of the dwarfs. The lack of very metal- 
rich and metal-poor stars can be explained by the fact that most of 
the stars are originated in the solar vicinity. Evolved stellar samples 
mostly consist of massive stars, which have shorter lifetime than 
the dwarfs, and therefore do not have enough time to migrate from 
further inner/outer disks ( [Wang & Zhao|2013][Minchev et al.|20r^ . 

Our present sample, as most of the giant star samples searched 
for planets, is affected hy B -V colour cut-off which excludes low- 
log g stars with high-[Fe/H] and high-log g stars with low metallic¬ 
ity. As discussed in |Mortier et al.[ ( [2013^ , this selection bias might 
be the reason of the absence of the correlation between occurrence 
of giant-planet planets and stellar metallicity. We suggest to use 
stars in a ”cut-rectangle“ in the log g - [Fe/H] diagram to overcome 
the aforementioned issue, if an unbiased sample is not available on 
hand. 

Although the current sample still contains only one star known 
to orbit a planetary companion (HD 11977 - [Setiawan et al.|2005^ 
most of the stars have already been periodically observed over the 
last years. Before a significant number of planets are detected, this 
sample can be used as a homogeneous comparison sample to study 
planet occurrence around giant stars. However, when exploring 
chemical peculiarities of planet-hosting giant stars, one should bear 


in mind the chemical properties of these evolved stars discussed in 
this paper (e.g., enhancement in Na, Al, etc.). 
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APPENDIX A: INTERDEPENDENCE OE STELLAR 
PARAMETERS AND THE MICROTURBULENCE 

The interdependence of the fundamental parameters are presented 
in Pig. |A1| The figure reveals several interesting correlations be¬ 
tween the parameters, for instance one can see that the metallicity 
correlates with surface gravity (see also Sect. 5) and also stars with 
higher Teff (above 5100 K) show higher metallicity. Microturbu- 
lent velocity correlates with the log g. The significance of the ob¬ 
served correlations are estimated following the method described 
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Table Al. The coefficients of the linear fits (y = a X X +b) of the relations 
between the stellar parameters, along with the correlation coefficient and 
the significance. The number of stars is 251. 


Elem 

a 

b 

R2 

z-score 

it -Teff 

0.120±0.065 

0.825±0.326 

0.013 

1.7 

logS-Teff 

0.847±0.089 

-1.356±0.441 

0.266 

7.9 

[Fe/H] - Teff 

0.280±0.069 

-1.472±0.343 

0.061 

3.9 

[Fe/H] - log g 

0.234±0.041 

-0.751±0.116 

0.116 

5.4 

it-log g 

-0.440± 0.029 

2.673±0.083 

0.476 

10.8 

it - [Fe/H] 

-0.154±0.057 

1.407± 0.010 

0.027 

2.6 


in [Figueira et al.| j2013| l and [Adibekyan et al.| ( |2013b and the pa¬ 
rameters of the linear relations are presented in Table |Al] We note 
that five stars classified as “outliers” in the ^t-log g, were excluded 
from the estimation of the significance of the correlations (see next 
section for details). 


Al The microturbulence relationship 

Sometimes, when the number of iron lines is not large enough, 
a correct determination of microturbulence becomes very difficult 
because of small EW interval of the Fel lines (e.g., |Mortier et ak] 
|2013| l. In these cases, one uses empirically obtained relations be¬ 
tween microturbulence and other stellar parameters. Several stud¬ 
ies have shown that for FGK dwarf stars, microturbulent velocity 


depends on log g and Tefi (e.g., Nissen|1981[[Allende Prieto et al. 

2004 

|Adibekyan et aL|2012||Tsantaki et al.|2013l Ramirez et al. 

2013 

1 . Takeda et al.|(2008|( has already suggested that the micro- 


turbulence correlates with the surface gravity, however the authors 
did not provide the analytic form of the relation. 


To find out the parameters the correlates with, we first ap¬ 
plied a linear fit for three pairs of data-sets: ^t-[Fe/H], ^t-logg. 
Then we evaluated the significance of the correlation, by 
using a bootstrap procedure as it was done in jFigueira et al.| ( (2013| l. 
As expected the strongest correlation is observed with log g (5.1a), 
~4(t in case of Tefi, and «1.8a for [Fe/H]. However, the fits can 
be affected by the presence of several outliers as can be seen in 
Fig. jAlj To remove the outliers we used the ^t-log g relation (since 
it shows the strongest correlation), by applying 2(T-clipping (two 
times of residual standard deviation). Then, after cleaning the data 
from outliers we again fitted the data and again evaluated the signif¬ 
icance of the relations. We found that microturbulence significantly 
correlated with the surface gravity (at about Her level), and with 
the metallicity but with less degree of significance. The five ouliers 
were responsible for the “strong” relation observed between and 

Teff- 

After this test, we decided to present the relation of microtur¬ 
bulence only with log g and [Fe/H], which has the following func¬ 
tional form: 

= 2.72(±0.08) - 0.457(±0.031) X logp 

-f0.072(±0.044) X [Fe/H] (Al) 

We note that this empirical relation is valid only for the range 
of stellar parameters that the stars in our sample cover. 



[a/Fe] 


Figure Cl. High-a and low-a separation histograms for the stars with 
metallicities < -0.3 dex (bottom), -0.3 Si [Fe/H] ^ 0.0 dex (middle), and 
[Fe/H] > 0.0 dex (top). The dotted lines are the separation curves between 
the thin and thick disks. 

APPENDIX B: [X/FE] DEPENDENCE ON STELLAR 
PARAMETERS 

In this section we present [X/Fe] vs. log g (Fig. [Blj l, [X/Fe] vs. 
microturbulence (Fig. |B2| (, and [X/Fe] vs. [Fe/H] (Fig. |B3^ plots 
derived from the “best” lines as discussed in the main text. 


APPENDIX C: SEPARATION OF THE GALACTIC DISKS 
BY q-ENHANCEMENT 

For the separation of Galactic stellar population by the chemical 
properties of the stars was done following the method presented 
in [Adibekyan et ar|p011| >. We first divided the sample into three 
metallicity bins: [Fe/H] < -0.3 dex, [Fe/H] > 0.0 dex, and stars 
in between. For the lowest and highest metallicity bins we easily 
identified the minima in the [a/Fe] histograms. For the interme¬ 
diate metallicity stars just plotting the [a/Fe] histogram will not 
reveal the minima, because the stars at these metallicities show a 
decrease of [a/Fe] with [Fe/H] (see Fig. [Cl[ l. Thus, we first de¬ 
trended the [a/Fe] by applying a liner fit and subtracting it. Then 
in the [a/Fe] histogram we identified the minima and by adding it 
to the previously applied liner fit we obtained the line which sep¬ 
arates the high- and low-a stars at -0.3 ^ [Fe/H] 5/ 0.0 dex. The 
separation lines for each metallicity bin presented in Fig.[Cl[ 
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[Fe/H] 


Figure Al. Interdependence of the stellar atmospheric parameters of the sample stars. The blue dotted lines depict the linear fits of the full data, and the sold 
lines are the fits after removing five “outliers”. The 2cr intervals of the linear fit of the relation are shown in blue dashed lines. The black croses 

indicate the five outliers. 
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Figure Bl. [X/Fe] vs. log g plots. Each element is identified in the upper right corner of the respective plot. The black dots represent the stars of the cun'ent 
sample. 
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Figure B2. [X/Fe] vs. microturbulence plots. Each element is identified in the upper right corner of the respective plot. The black dots represent the stars of 
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Figure B3. [X/Fe] vs. [Fe/H] plots derived from the “best” lines. For Nal, CrII, and Mnl there was no “best” line(s) foud. The black dots represent the stars 
of the sample and the gray small dots represent stars from |Adibekyan et al. 1^201^ with Teff = Te ±500 K. The red circle and blue square show the average 
[X/Fe] value of stars with [Fe/H] = 0.0±0.1 dex. Each element is identified in the upper right corner of the respective plot. 


















































